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Abstract  Ocean  surface  waves  (sea  and  swell)  are  generated  by  winds  blowing  over  a  distance  (fetch)  for 
a  duration  of  time.  In  the  Arctic  Ocean,  fetch  varies  seasonally  from  essentially  zero  in  winter  to  hundreds 
of  kilometers  in  recent  summers.  Using  in  situ  observations  of  waves  in  the  central  Beaufort  Sea,  combined 
with  a  numerical  wave  model  and  satellite  sea  ice  observations,  we  show  that  wave  energy  scales  with  fetch 
throughout  the  seasonal  ice  cycle.  Furthermore,  we  show  that  the  increased  open  water  of  201 2  allowed 
waves  to  develop  beyond  pure  wind  seas  and  evolve  into  swells.  The  swells  remain  tied  to  the  available 
fetch,  however,  because  fetch  is  a  proxy  for  the  basin  size  in  which  the  wave  evolution  occurs.  Thus,  both  sea 
and  swell  depend  on  the  open  water  fetch  in  the  Arctic,  because  the  swell  is  regionally  driven.  This  suggests 
that  further  reductions  in  seasonal  ice  cover  in  the  future  will  result  in  larger  waves,  which  in  turn  provide  a 
mechanism  to  break  up  sea  ice  and  accelerate  ice  retreat. 


1.  Introduction 


In  the  open  ocean,  surface  waves  (sea  and  swell)  are  a  central  feature  of  the  air-water-ice  interface,  where 
significant  exchanges  of  momentum,  heat,  and  gases  occur  [Steele  et  at.,  1 989;  Melville,  1 996].  At  the  coasts, 
surface  waves  can  force  circulation  and  cause  erosion  [ Battjes ,  1988].  Surface  waves  also  can  be  a  haz¬ 
ard  to  navigation  and  operations  at  sea.  In  the  emerging  Arctic  Ocean,  all  of  these  wave  processes  are 
becoming  significant  to  the  natural  and  human  systems  [Stephenson  et  ai,  2011;  Showstack,  2013].  The 
interaction  of  waves  and  ice  is  particularly  complex,  because  ice  can  suppress  waves  by  scattering  and 
dissipating  wave  energy  [Squire  et  ai,  1 995;  Squire,  2007],  while  the  waves  simultaneously  break  up  the 
ice  [Dozaki  etal.,  1999;  Marko,  2003].  Ice  also  has  an  indirect  effect  on  surface  waves,  by  setting  the  fetch 
available  for  wave  generation  and  the  distance  available  for  wave  evolution.  This  is  the  effect  we  study  in 
this  paper. 

In  the  absence  of  ice,  the  dependence  of  ocean  surface  waves  on  fetch  is  a  well-studied  problem,  and 
fetch  has  been  found  to  be  the  limiting  variable  for  wave  generation  in  coastal  regions  and  marginal  seas 
[Hasselmann  etal.,  1973],  Several  empirical  fits  describe  the  evolution  of  waves  along  a  fetch  x  acted  on  by  a 
uniform  wind  U  using  the  nondimensional  scaling  parameters  of  wave  energy  and  fetch  [Young,  1999]: 
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where  H  is  significant  wave  height  and  g  is  gravity.  Storm  duration  can  be  similarly  scaled  as  T  =  jj,  and 
duration  is  often  the  limiting  variable  for  wave  growth  in  the  open  ocean  [Young,  1999].  In  practice,  uni¬ 
form  winds  are  rarely  observed,  yet  these  existing  empirical  fetch  laws  often  remain  in  agreement  with 
observations  [Schwendeman  et  ai,  2014],  For  a  pure  wind  sea  (i.e.,  short  waves  forced  by  local  winds),  the 
Pierson-Moskowitz  limit  asserts  a  maximum  nondimensional  energy  £  =  3.64  x  10-3  that  occurs  for  a 
saturated  sea  state  [Pierson  and  Moskowitz,  1964].  This  limit  remains  debated  and  undervalidated,  in  part 
because  most  natural  wave  conditions  include  both  sea  and  swell  (i.e.,  long  waves  resulting  from  nonlocal 
winds)  [Alves  et  ai,  2003].  Swell  waves  accumulate  the  energy  in  the  wind  sea  and,  given  enough  space  and 
time,  often  exceed  this  limit. 


Flere  we  intentionally  include  conditions  with  both  sea  and  swell  and  examine  the  large-scale  effects  of 
variable  open  water  distances  in  the  Arctic  Ocean.  This  is  in  contrast  to  the  strict  definition  of  a  fetch  depen¬ 
dence,  in  which  a  pure  wind  sea  achieves  a  local  equilibrium  between  generation  by  wind  forcing  and 
dissipation  by  wave  breaking  [Fontaine,  2012],  Rather,  we  use  the  distance  scale  X  to  describe  the  overall 
space  available  for  waves  to  be  generated  as  a  wind  sea  and  then  evolve  to  become  swell. 
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In  the  modern  Beaufort  Sea,  the  open 
water  distances  change  dramatically 
throughout  the  summer  season,  from 
essentially  zero  in  April  (sea  ice  cover 
maximum)  to  well  over  1000  km  in 
September  (sea  ice  cover  minimum).  In 
recent  years,  the  seasonal  ice  retreat  has 
expanded  dramatically  [Stroeve  et  al., 
2005,  2008;  Simmonds  and  Keay,  2009], 
leaving  much  of  the  Beaufort  Sea  ice 
free  at  the  end  of  summer.  For  exam¬ 
ple,  on  1 8  September  201 2  strong  winds 
(up  to  1 8  m  s_1 )  generated  5  m  wave 
heights  in  the  ice-free  central  Beaufort 
Sea  (Figure  1).  These  conditions  were 
measured  in  situ  using  a  mooring  near 
the  middle  of  the  basin  and  modeled 
using  a  spectral  wave  action  balance 
code  (WAVEWATCH  III  ®).  There  are  few 
previous  studies  of  waves  in  the  Arctic 
Ocean;  they  have  focused  mostly  near 
the  coasts  or  on  local  effects  [Francis  et 
al.,  2011;  Squire  et  al.,  2009;  Bogucki  et  al., 
2013].  Of  the  previous  studies,  Francis  et 
al.  [201 1]  in  particular  showed  a  strong 
trend  of  increasing  surface  wave  heights 
in  the  Arctic.  The  observations  reported 
here  are  the  only  known  wave  mea¬ 
surements  in  the  central  Beaufort  Sea, 
because  until  recently  the  region  remained  ice  covered  throughout  the  summer  and  there  were  no  waves 
to  measure. 

2.  Methods 

Wave  observations  were  collected  using  a  600  kFIz  Nortek  Acoustic  Wave  and  Current  (AWAC)  sensor 
deployed  on  a  subsurface  mooring  at  75°  N,  150°W  as  part  of  the  Beaufort  Gyre  Exploration  Project.  The 
AWAC  was  mounted  at  the  top  (50  m  below  the  surface)  of  the  3000  m  mooring,  looking  up  to  the  surface. 
The  AWAC  uses  sonar  pings  to  determine  the  distance  from  the  instrument  to  the  sea  surface,  and  this  dis¬ 
tance  is  recorded  at  1  Flzfor  1024  s  at  the  beginning  of  each  hour.  These  records  are  processed  using  a  fast 
Fourier  transform  to  determine  the  wave  energy  spectra,  as  well  as  the  significant  wave  height  (defined  as  4 
times  the  square  root  of  the  variance  in  sea  surface  elevation  at  frequencies  0.05-0.6  Flz  [Me/,  1989]),  or  the 
ice  draft  when  ice  is  present  [Magnell  et  al.,  201 0]. 

Wave  hindcasts  were  made  using  the  WAVEWATCH  III®  numerical  wave  model  [Tolman,  2009],  forced  with 
the  U.S.  Navy's  Fleet  Numerical  Meteorology  and  Oceanography  Center  (FNMOC)  wind  fields  and  ice  prod¬ 
ucts  (derived  from  satellite  passive  microwave  radiometers).  Ice  presence  or  absence  was  determined  using 
a  25%  threshold  [Tolman,  2003]. 

Fetch  estimates  were  made  using  the  NSIDC  [National  Snow  and  Ice  Data  Center,  201 2]  4  km  daily  ice  prod¬ 
uct  and  calculating  the  distance  along  the  direction  of  the  FNMOC  winds  from  the  mooring  location  to 
an  ice  or  land  boundary.  During  the  maximum  open  water  period  (September),  the  mooring  is  located 
such  that  the  fetch  distance  exceeds  1 000  km  in  all  directions  and  the  fetch  calculation  is  a  proxy  for  the 
basin  scale. 

Storm  duration  estimates  were  made  as  the  cumulative  hours  during  which  FNMOC  wind  speed  exceeded 
5  m/s,  with  a  reset  to  zero  whenever  the  wind  dropped  below  5  m/s. 


0  0.5  1  1.5  2  2.5  3  3.5  4  4.5  5 

Wave  height  [m] 


Figure  1 .  Example  wave  model  hindcast  during  September  201 2  storm. 
The  map  is  centered  on  the  North  Pole,  and  the  mooring  location  is 
indicated  by  the  black  circle  north  of  Alaska.  The  color  scale  indicates 
significant  wave  height  from  0  to  5  m. 
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Figure  2.  Time  series  of  Arctic  Ocean  wave  conditions.  Hourly  values  for  (a)  significant  wave  height,  (b)  wind  speed  at 
10  m  reference  height,  and  (c)  fetch  along  the  direction  of  the  wind,  from  the  mooring  location  to  ice  or  land.  In 
Figure  2a,  symbols  show  the  in  situ  observations,  and  the  solid  curve  shows  the  model  hindcast  results.  Gray  is  used 
to  indicate  the  presence  of  partial  ice  cover  (as  determined  from  the  mooring  data).  In  Figure  2c,  filled  symbols  are  for 
wave  conditions  that  are  pure  wind  seas  and  open  symbols  are  for  wave  conditions  that  include  swells  (as  determined 
by  nondimensional  wave  age,  ^). 


3.  Results 

There  is  excellent  agreement  in  wave  height  between  the  observations  for  most  of  the  summer  2012  record 
and  the  model  (Figure  2),  but  there  is  a  notable  bias  in  mid-August  (i.e.,  immediately  following  the  mooring 
deployment  on  1 2  August)  when  the  model  overestimates  the  observations  by  as  much  as  a  factor  of  2.  This 
corresponds  to  the  seasonal  retreat  and  advance  of  sea  ice,  during  which  the  NSIDC  products  indicate  that 
partial  cover  occurs  at  the  mooring  location.  The  wave  model  includes  methods  for  treating  the  effect  of 
ice  (using  estimates  of  ice  concentration  as  input),  but  it  is  believed  that  the  ice  products  employed  here, 
derived  from  satellite  passive  microwave  radiometer,  underrepresent  the  ice.  Clearly,  the  local  effects  of  ice 
on  the  waves  at  a  given  location  are  significant  and  warrant  further  investigation.  Here,  however,  we  focus 
on  the  regional  effects  of  sea  ice  variations  and  examine  the  middle  portion  of  the  wave  record,  when  the 
mooring  location  is  ice  free  but  ice  cover  farther  north  of  the  mooring  location  continues  to  control  the 
effective  size  of  the  basin. 

There  were  two  notable  wave  events  in  2012:  a  storm  early  August  and  another  in  mid-September.  The 
first  storm  occurred  before  the  in  situ  observations  began  (and  when  satellite  observations  suggest  that 
ice  cover  was  still  significant  at  the  mooring  location);  hence,  only  model  results  are  shown  for  that  storm. 
Although  the  first  storm  was  significant  in  forcing  ice  retreat  [Parkinson  and  Comiso,  2013],  the  waves  asso¬ 
ciated  with  that  storm  were  not  modeled  to  be  as  large  as  the  September  storm.  This  is  likely  because  the 
fetch  available  for  wave  growth  was  much  smaller  in  August  (when  the  ice  edge  was  at  approximately  72°N) 
than  in  September  (when  the  ice  edge  was  at  approximately  80°N). 

The  overall  dependence  on  open  water  distance  is  shown  in  Figure  3,  where  nondimensional  wave  energy 
scales  with  nondimensional  distance  over  3  orders  of  magnitude.  (The  points  from  the  September  storm, 
for  example,  lie  near  the  regression  line  and  span  104  to  1 06.)  The  coefficient  of  determination  is  R 2  =  0.81, 
which  is  significant  beyond  the  95%  level  and  is  notably  higher  than  the  direct  correlation  of  wind  speed 
and  wave  height  (R2  =  0.67,  not  shown).  In  contrast,  there  is  no  significant  relation  between  nondimen¬ 
sional  wave  energy  and  nondimensional  storm  duration  (R2  =  0.001,  not  shown). 

Scaling  analysis  such  as  this  is  often  prone  to  spurious  correlation,  because  the  wind  speed  U  is  used  in  both 
variables.  Here  the  lack  of  correlation  between  wave  energy  and  duration,  which  also  each  have  U  in  their 
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scaling,  confirms  that  the  fetch  correla¬ 
tion  is  not  spurious.  Furthermore,  there  is 
no  systematic  dependence  or  correlation 
with  wind  direction,  which  is  distributed 
nearly  uniformly  over  the  record  in 
the  range  0-15  m  s-1.  Other  relation¬ 
ships,  however,  are  apparent  in  the  data. 
Dimensionless  peak  period  has  a  strong 
positive  correlation  with  dimensionless 
fetch  ( R 2  =  0.7)  and  a  dependence  A*1/2. 
Dimensionless  energy  has  a  strong  cor¬ 
relation  ( R 2  =  0.5)  with  wave  age  and  a 
dependence  (c/L/)4. 
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Figure  3.  Scaling  of  waves  in  the  Arctic  Ocean,  using  nondimensional 
wave  energy  versus  nondimensional  fetch.  There  are  1880  points  (each 
corresponding  to  an  hourly  observation  from  the  in  situ  mooring  at 
75°N,  150°W).  Symbols  as  in  Figure  2.  The  dashed  line  is  a  regression 
with  a  logarithmic  slope  of  1.6.  The  Pierson-Moskowitz  limit  for  pure 
wind  seas  is  shown  at  S  =  3.64  x  10~3. 


The  empirical  scaling  (Figure  3)  obscures 
many  of  the  details  of  the  wave  action 
balance  and  wave-ice  interactions. 

In  contrast  to  pure  wind  sea  studies 
[Dobson  et  at.,  1 989;  Stiassnie,  201 2], 
the  scaling  includes  swell  (e.g.,  domi¬ 
nant  wave  periods  of  up  to  1 0  s  during 
the  September  storm),  and  the  values 
greatly  exceed  the  Pierson-Moskowitz 
limit.  The  swell  is  not  driven  by  the  local 

wind  and  thus  is  an  external  source  of  wave  energy  in  the  scaling.  To  isolate  points  with  swell,  the  points 
in  Figure  3  are  separated  by  nondimensional  wave  age,  which  is  the  ratio  of  the  dominant  wave  celer¬ 
ity  to  the  wind  speed.  Pure  wind  seas  have  a  wave  age  less  than  one,  indicating  that  the  wind  is  driving 
the  waves,  and  these  points  cluster  largely  below  the  Pierson-Moskowitz  limit.  Swells  have  a  wave  age 
greater  than  one,  indicating  that  the  waves  are  outrunning  the  wind,  and  these  points  extend  well  above 
the  Pierson-Moskowitz  limit.  The  seas  have  the  expected  dependence  of  X1 ,  whereas  the  swells  have  a 
dependence  closer  to  A*1  7.  The  overall  fit  is  A’1-6. 


The  scaling  is  successful  across  both  categories,  because  these  waves  represent  a  continuum  of  sea  and 
swell.  The  synoptic  wind  fields  that  drive  the  wind  sea  are  large  enough  in  scale  to  be  relevant  to  the  whole 
basin  as  the  wind  sea  evolves  into  swell.  This  is  similar  to  the  results  of  Young  [2006],  in  which  both  sea  and 
swell  follow  a  fetch  law  in  Flurricanes  (his  Figure  11)  because  under  extreme  forcing  the  swell  remains  part 
of  the  wind  and  wave  system.  Flere  the  distance  of  open  water  is  relevant  to  the  nonlinear  interactions  that 
control  the  development  of  swells  (because  the  nonlinearity  is  weak  and  requires  large  propagation  dis¬ 
tances)  [Hasselmann,  1962].  In  the  limit  of  pure  swell,  the  local  wind  is  not  important,  but  the  scaling  of  wave 
height  and  basin  size  is  still  valid  because  space  is  required  to  build  the  swell.  For  swells,  the  length  scale  X 
can  be  interpreted  as  the  basin  scale  (although  the  actual  fetch  is  still  used),  because  a  mooring  in  the  center 
of  the  basin  will  have  a  similar  value  (regardless  of  wind  direction)  that  represents  the  characteristic  distance 
to  an  ice  boundary. 


4.  Conclusion 

The  Arctic  Ocean  is  unique  in  annual  variation  of  length  scale,  with  the  maxima  increasing  dramatically 
in  recent  years.  Although  the  dynamics  are  more  general  than  the  conventional  fetch  dependence  of 
ocean  waves,  the  result  is  simple:  open  water  distances  are  the  controlling  variable  for  wave  heights  in  the 
Beaufort  Sea  (and  likely  the  rest  of  the  Arctic  Ocean).  Future  scenarios  for  reduced  seasonal  sea  ice  cover  in 
the  Arctic  suggest  that  larger  waves  are  to  be  expected  and  that  swells  will  be  more  common.  Swells  carry 
more  energy  and  have  longer  attenuation  scales  within  ice  [Squire  et  at.,  1 995;  Squire,  2007]  and  thus  will 
be  more  effective  at  breaking  up  the  remaining  ice.  It  is  possible  that  the  increased  wave  activity  will  be  the 
feedback  mechanism  which  drives  the  Arctic  system  toward  an  ice-free  summer.  This  would  be  a  remark¬ 
able  departure  from  historical  conditions  in  the  Arctic,  with  potentially  wide-ranging  implications  for  the 
air-water-ice  system  and  the  humans  attempting  to  operate  there. 
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The  increasing  wave  climate  will  also  have  implications  for  the  coasts,  which  are  already  eroding  rapidly  dur¬ 
ing  summer  months  as  a  result  of  climate  change  and  subsequent  loss  of  permafrost  [Overeem  et  a!.,  2011], 
Using  the  conventional  approach  that  shoreline  sediment  transport  depends  primarily  on  incident  wave 
energy  [Bailard,  1981],  the  expectation  then  is  for  coastal  erosion  to  increase  with  a  similar  A’1  6  dependence 
on  the  open  water  fetch  of  the  Arctic  Ocean.  Although  the  actual  wave  effects  will  of  course  be  site  specific 
and  complex,  our  scaling  is  a  starting  point  in  understanding  the  rapidly  changing  wave  climate  in  the  Arctic 
Ocean  and  the  likely  expanding  future  role  of  waves  in  the  Arctic  system. 
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